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Abstract: We study the kinematic Sunyaev-Zel'dovich (kSZ) effect in a Lemitre-Tolman- 
Bondi (LTB) universe model whose distance-redshift relation agrees with that of the con- 
cordance ACDM model at redshifts z <2. This LTB universe model has a void with size 
comparable to the Hubble horizon scale. We first determine the decoupling epoch in this 
LTB universe model by an approximate analytical condition under a few simplified assump- 
tions on the physical quantities at that epoch. Then we calculate the cosmic microwave 
background (CMB) anisotropy observed in the rest frame of clusters of galaxies which are 
assumed to be at rest in the spatial comoving coordinates of the LTB universe model. We 
find that the obtained temperature anisotropics are dominated by dipole, although there 
may exist higher multi-poles in general. We may interpret this dipole anisotropy as the 
drift velocity of a cluster of galaxies relative to the CMB rest frame. Hence it gives rise to 
the kSZ effect. We calculate this effect and compare it with observational data. We find 
that if we assume the conventional adiabatic perturbation scenario at the time of decou- 
pling, the drift velocity of clusters of galaxies becomes unacceptably large. Conversely, this 
observational constraint may be relaxed by introducing a non-adiabatic (i.e., primordially 
isocurvature) component of inhomogeneities at the time of decoupling. However, our result 
indicates that the necessary isocurvature perturbation amplitude is very large. 



Contents 



|l]. Introduction ^ 

|2|. Decoupling Epoch in the LTB Universe § 

|2.1| Background universe model ^ 

Decoupling epoch in the LTB universe S 



|3|. Kinematic Sunyaev-Zel'dovich effect ^ 
|3.1| Case of homogeneous decoupling hyper surf ace 10 



3.2 Case of inhomogeneous decoupling hypersurface 11 



HL Summary and Discussion 14 



1. Introduction 

Copernican principle, i.e., the assumption that we are not at a special place in the universe, 
plays a very crucial role in the modern physical cosmology. If we accept this principle, 
the observed highly isotropic distribution of the cosmic microwave background (CMB), 
except for its dipolar anisotropy, implies the homogeneity and isotropy of our universe. 
The standard cosmological model based on this homogeneity and isotropy and on general 
relativity, namely the ACDM model, naturally explains various important observational 
facts, supporting strongly the validity of the Copernican principle, though the homogeneity 
of our universe has not yet been directly tested. 

Within the framework of homogeneous and isotropic cosmological models, the distance- 
redshift relation of type la supernovae and the CMB data by WMAP lead us to a conclusion 
that there is so-called dark energy as a major component of the universe. The present 
observational data is consistent with the dark energy being a positive cosmological constant. 
However, the energy density of the dark energy is about 10^^^*^ times smaller than the 
Planck energy density. This number is regarded as one of the most unnatural numbers of 
the universe. Hence, recently, many authors have been discussing the possibility to explain 
observations without dark energy by using anti-Copernican universe models in which we are 
assumed to be at a very special position in the universe. Usually, in anti-Copernican models, 
we are assumed to be located at at the center of a spherically symmetric inhomogeneous 
universe |,|,|, |,§. 

One of the important attempts in this direction is to explain the type la supernovae 
observation using spherically symmetric inhomogeneities without dark energy. There exits 
an exact solution of the Einstein equations known as the Lemaitre-Tolman-Bondi (LTB) 
solution which describes a system of spherically symmetric dust. Actually, it is possible to 
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construct an LTB universe model whose distance-redshift relation observed at the symme- 
try center agrees with that of the concordance ACDM model ^ 0, ^, |9|, |lO|. Such an 
LTB universe model has a very large void structure comparable to the cosmological horizon 
scale. This result is very important because it implies that, in order to know the precise 
information about the dark energy, we have to know how to distinguish the effect of the 
dark energy from that of inhomogeneity. For this purpose, in this paper, we also study the 
LTB cosmology with an observer like us at its center. 

CMB observations in LTB universe models are also often discussed 0, H, 

1^, H, m, i0|, H, H, H, H, H, H, m. in many of these works, it is 



assumed that the universe is homogeneous in the spacelike asymptotic region from which 
the CMB photons come, hence the CMB photon distribution is assumed to be the same as 
that in the homogeneous and isotropic universe models at the last scattering surface (LSS). 

One of the most stringent constraints on the LTB universe model with a large void 
come from the observation of the kinematic Sunyaev-Zel'dovich (kSZ) effect. In the case 
of a homogeneous and isotropic universe, if a cluster of galaxies has a finite drift velocity 
relative to the CMB comoving frame, it causes a dipole anisotropy in the CMB temperature 
in the cluster rest frame. Then Compton scattering of the CMB photons inside a cluster of 
galaxies causes a distortion in the CMB spectrum in the direction of the cluster of galaxies. 
Therefore observation of distortions in the CMB spectrum constrains the drift velocity of 
clusters of galaxies. 

In the case of LTB universe models, even if clusters of galaxies are comoving with the 
non-relativistic matter component in the universe, or in other words, even if they are at 
rest in the comoving spatial coordinates, an anisotropy of the CMB temperature will be 
observed at an off-center cluster due to the inhomogeneity of the universe. In contrast 
to a homogeneous and isotropic universe, this effect may not necessarily be of kinematic 
origin, since higher multi-pole components of the CMB anisotropy may exist in the cluster 
rest frame. This is because an LTB universe model is neither homogeneous nor isotropic 
at off-center clusters of galaxies. In this sense, the terminology 'the kSZ effect' is not so 
appropriate for this effect in the LTB universe model. However, as will be shown later, in 
our LTB universe model, the anisotropy observed in the rest frame of a cluster of galaxies is 
found to be dominated by a dipole component. Since a dipole anisotropy may be regarded 
as the effect of the drift velocity of a cluster of galaxies, we may call this effect the kSZ 
effect. 

In Ref. [0], Garcia-Bellido and Haugboelle reported that current observations of only 
nine clusters with large error bars would already rule out an LTB universe model with 
a void of a radius greater than ~ 1.5Gpc. Their LTB universe model also has a large 
spherical void structure, but it is homogeneous and isotropic in the spatially asymptotic 
region. Thus, the CMB at and near the last scattering surface is estimated by the same 
procedure as in the case of a homogeneous and isotropic universe model. Differently from 
standard cosmology, however, in the LTB universe models, it might be possible to introduce 
radial inhomogeneities in the CMB temperature. 

In this paper, we analyze the kSZ effect in the LTB universe model whose distance- 
redshift relation observed at the center agrees with the concordance ACDM model, and 
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then we show that the observational constraints from the kSZ effect have ruled out LTB 
universe models based on the adiabatic perturbation scenario, but not so if we allow the 
existence of non-adiabatic inhomogeneities at the decoupling epoch. 

This paper is organized as follows. In §|2|, we introduce the LTB solution and show how 
to fix the decoupling epoch in an LTB universe model. Under some reasonable assumptions 
on the physical quantities at the decoupling epoch, we calculate the dipole anisotropy of 
the CMB temperature observed at each cluster of galaxies in the adiabatic scenario in 
which there existed only adiabatic curvature perturbations in the early universe in § |3.1| . 
We find that the dipole anisotropy is too large to be consistent with observational data. 



Then in § |3.2| , we assume the existence of isocurvature inhomogeneities at the decoupling 
epoch, and show that it may save LTB universe models. Finally, §^ is devoted to summary 
and discussion. 

2. Decoupling Epoch in the LTB Universe 
2.1 Background universe model 

In Ref. Yoo, Kai and Nakao numerically constructed an LTB universe model whose 
distance-redshift relation agrees with that of the concordance ACDM model in the whole 
redshift domain and which is homogeneous at the early stage of the universe. Recently 
we modified this LTB universe model at redshifts z > 2 so that the peak positions in the 
CMB temperature anisotropy are consistent with the observation [|3l|] (hereafter referred 
to as Paper I). In this paper, we use this LTB model as the background universe. 
The metric of the LTB solution is given by 

ds^ = -c^de + ^2 ^ ^2(^^^)(^^2 ^ sin2^d(^2)^ (2.1) 

1 — k[r)r^ 

where k{r) is an arbitrary function of the radial coordinate r. The matter is dust whose 
stress-energy tensor is given by 

T^"' = pu'^u", (2.2) 

where p = p{t, r) is the mass density, and is the four-velocity of the fluid element. The 
coordinate system in Eq. ( |2.1| ) is chosen in such a way that = (1,0,0,0). 
The area radius R{t, r) satisfies one of the Einstein equations, 

^^V = ^-.W, (2.3) 



dt J R 

where M{r) is an arbitrary function related to the mass density p by 



. , 1 dMir) , , 



Following 1 32], we write the solution of Eq. ( |2.3| ) in the form, 

R{t, r) = (6GM(r))i/3(t - tB{r)fl^ S{x), (2.5) 
= ^k{Ty (2.6) 



V 6GM(r) 
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where tB{r) is an arbitrary function which determines the big bang time, and S{x) is a 
function defined imphcitly as 



S{x) 



cosh ^/—r] — 1 



- (sinh y/—r] 



17^)2/3 



6-*^/^ (sinh ^/—r] — yj—r])'^/^ 

1 — cos 

6^/3 — sin -y/r/)^/^ ' 



62/3 

- sin 



62/3 



for X < . 
for X > . 



(2.7) 



and 5(0) = (3/4)-*^/^. The function S{x) is analytic for x < (tt/S)^/^. Some characteristics 
of the function S{x) are given in ||3^ and Q. 

As shown in the above, the LTB solution has three arbitrary functions, k{r), M{r) 
and tB{r)- One of them is a gauge degree of freedom for rescaling of the radial coordinate 
r. We fix this by setting 

4 



M(r) 



-npor 



2.8 



where po is the energy density at the center at present po = p{to, 0). As in the case of the 
homogeneous and isotropic universe, the present Hubble parameter Hq is related to po as 



Hi + k{Oy = -ttGpo. 
As in |8|, we assume the simultaneous big bang, i.e.. 



tB{r) = 0. 



(2.9) 



(2.10) 



For notational simplicity, we introduce dimensionless quantities, 

Hor ~ k{r)c^ 
r:=—-, k{r):=—^. 

The LTB universe model proposed in Paper I is specified by the following curvature 
function: 



k{^ = A;fit(F) X f{r;A), 



where 



k&tir) 



0.545745 



2.22881 



0.211472 + V0.026176 + f (0.807782 + V0.026176 + f) 
1 for r < 2 

^ , 16A(F-2)3(323-123F+12?^) „ n /n 

1 H ^ 3125 2 < r < 9/2 

l + A for 9/2 < r 



-\2 ' 



(2.11) 

(2.12) 
(2.13) 



with A = —1.069. The distance-redshift relation of this LTB solution agrees with that of 
the concordance ACDM model at redshifts z 2. 
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2.2 Decoupling epoch in the LTB universe 

The decoupHng between photons and baryons occurs in an inhomogeneous universe just 
as in the case of a homogeneous universe. That is, it occurs when the mean free path of 
photons becomes effectively infinite due to almost complete recombination of electrons to 
protons. 

Since there is no radiation component in our LTB model, we cannot treat the decou- 
pling in a rigorous manner. However, as in the concordance model, we expect the energy 
density of the radiation to be only a small fraction of the total density at decoupling, hence 
its effect on the spacetime geometry is small, if not negligible. In fact, the radiation energy 
density estimated in our LTB model turns out to be about 20% of the total energy density. 
This means that treating the radiation as a test field in our model is consistent to a first 
approximation. 

Another approximation we adopt is the instantaneous decoupling. Namely, we assume 
decoupling to occur on a single spacelike hypersurface. Since our LTB universe model 
is inhomogeneous but spherically symmetric, it is natural to assume that the decoupling 
hypersurface is also inhomogeneous but spherically symmetric. Thus it is specified by the 
form, 

t = tB{r) . (2.14) 

The cross section of this hypersurface with the past directed null cone from the observer 
at the center constitutes the LSS of CMB photons (see Fig. |l|). The LSS is a spacelike 



Observer 




(sphere) 



Figure 1: Schematic figure for the hypersurface of the decoupling epoch in an LTB universe. 
2-dimensional sphere by the assumed symmetry. We use the subscript * to express phys- 
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ical quantities on the LSS. In our approximation, ignoring secondary effects, the CMB 
anisotropy is essentially determined by the distribution of photons on the LSS. 

The geodesic equations to determine the past light cone from the observer at the center 
are written in the form, 

, , dt OrR , , 



where the past directed radial null geodesies have been parametrized by the cosmological 
redshift z. We denote the solution of the above equations by 

t = tic(z) , r = rx^{z). (2.17) 

Now to discuss the decoupling condition, for simplicity, we consider a universe consists 
of cold dark matter, protons and electrons, and neutral hydrogen atoms. In particular, we 
neglect helium. Since the contributions of helium and the other components are not large, 
this simplification should not lead to a serious error in our analysis. We also assume that, 
until the decoupling time, photons, electrons and protons are in thermal equilibrium. The 
energy density of the electrons and photons is negligible, and hence the constituents of our 
LTB model are cold dark matter and baryons, the latter of which consist of protons and 
hydrogen atoms. Thus the baryon number density nb is equal to the total number density 
of protons and hydrogen atoms, and the electron number density ne is equal to the proton 
number density. 

In homogeneous and isotropic cosmology, the decoupling time is well determined by 
Gamow's criterion, 

F = r, (2.18) 

where H is the Hubble parameter and F is the rate of collisions of a photon with electrons. 
Using the Thomson scattering cross section ctt and the electron number density rie, T is 
written as 

r = cnef7T. (2.19) 

In our LTB model, we adopt this criterion, with the identification of the "Hubble parameter 
F" with 

^^ = ^P. (2.20) 

Since, by virtue of Eq. ( |2.10| ), our LTB universe model is almost identical to the Einstein- 
de Sitter universe near the decoupling time. Hence the above definition of H is accurate 
enough for our purpose. 

In order to estimate the electron number density, for simplicity, we use Saha's equation 
assuming thermal equilibrium. We note that our assumption of thermal equilibrium at the 
decoupling epoch is not appropriate in reality, because the decoupling process is controlled 



by a non-thermal physical process [^3[. However, as mentioned in |28|, the error in the 
temperature estimate turns out to be 2 - 3%. Hence our treatment of the decoupling is 
good enough for our qualitative discussions about the kSZ effect. 
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Let us consider the ionization rate := n^/n\^. In thermal equihbrium, the ionization 
rate X^, satisfies Saha's equation, 

4^/2C(3) (k^T\"^ /13^9eV\ 

where ^(ic) is the zeta function, and T, r/, k-Q and me are the temperature, the baryon- 
to-photon ratio, the Boltzmann constant and the electron mass, respectively. Since the 
ionization rate at decoupling drops down to X^. ~ 10~^, we may approximate the above 
equation by 

Xl . f^V^xp f-^) . (2.22) 

' 4^/2C(3)?? V"^ec2y ^ \ k^T ) ^ ' 

Using this equation, Gamow's criterion ( |2.18D is rewritten in the form, 

32v^C(3) Gp fkBToV^^ fny/' f l3.59eV \ 

where n^o and Tq ~ 2.725K are the present photon number density and observed CMB 
temperature, respectively. 

Since we assume thermal equilibrium of electrons, protons and photons until the decou- 
pling time t = t]:,{r), the physical state of the decoupling hypersurface, which is spherically 
symmetric, is determined by the distributions of the temperature T = To^r), the baryon- 
to-photon ratio r] = r/D(?'), and the matter energy density p = p{tD{r),r). For convenience, 
in place of p{ti:){r),r), we introduce the following quantity: 

Mr):=^MMf_Z^y. (2.24) 



PO V^D(r) 

The quantity aj){r) is proportional to the ratio of the matter density and the photon 
number density. Then, from Eq. (2.23), we obtain 

rj^(r) = 32V2^C(3) Gpo ( kBToV^^^^^^ ( J±_Y^\^p /13^59eV\ 25) 

We note that once Tb(r) and aD('") are given, the hypersurface t = t]:){r) can be obtained 
from Eq. p3^ ). 

On the LSS, we must have 

Tpinciz^)) T* 

7f = — = l + z*. (2.26) 

Then, if we regard ?b(r), rjDir) and aB{r) as mutually independent functions, we have one 
functional condition ( p. 25]) and one boundary condition at z = z^, to constrain these three 
functions. But these are not enough to determine the decoupling hypersurface, t = tnir), 
through Eq. ( |2.24 ). We would need two more functional conditions. However, to know the 
location of the LSS for the observer at the center, we only need a single condition on the 
LSS in addition to Eqs. ( 2.25 ) and ( |2.26| ). For simplicity, we impose 

r/, = 6.2 X 10"^° (2.27) 
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as in Paper I 

We numerically solve the radial null geodesic equations ( 2.15| ) and ( 2.16| ). At each 
redshift z, we define T and a by 

T = (1 + z)To, (2.28) 
a = P^'^^['^^'^^(^!^\ (2.29) 



During the numerical integration, we check at each redshift z whether Eq. ( p. 25 ) is satisfied 



by Td = T, qd = a and r/D = ??*• If Eq. ( 2.25 ) is satisfied, we stop integrating Eqs. ( 2.15| ) 



and (|2.16| ), and identify z, T and a at this moment with z*, and a*, respectively. For 
our LTB universe model, we have obtained 

^ = l + z, - 1130, (2.30) 
Jo 

a=,~ 5.184, (2.31) 
3. Kinematic Sunyaev-ZeI'dovich effect 

We consider the kSZ effect caused by the hot gases bound within clusters of galaxies which 
are assumed to be at rest in the spatial coordinates of the LTB universe model. Let 
us consider a cluster of galaxies at redshift z = Zc- In order to know how the CMB is 
observed at the position of this cluster of galaxies, we numerically integrate pastward the 
null geodesic equations from this cluster to every direction (see Fig. ^). Without loss of 
generality, we may focus on null geodesies in the plane ip = 0. Then the non-radial null 
geodesic equations are given by 

1 — Kr^ 

.2\ 



r" = -2 [dtilndrR)] t'r' + ^ [a^ln(l - kr"^) - 2dr{\ndrR)] r'^ + 1?'^(3.2) 
d" = 2{dt In R) t'd' - 2{dr In R) r'l?', (3.3) 



where the prime denotes differentiation with an affine parameter. From the null condition, 
we have 

r" = 1^ (t" - R' . (3.4) 

(drR)^ V ; 

The initial conditions (tj, n, t', r[, i?J) to integrate the null geodesic equations are given 

by 

ti = tc := tic{zc), (3.5) 

n = rc := ric{zc), (3.6) 

= 0, (3.7) 

t: = -(l + ze), (3.8) 

K = ^5^^(l + ^c)cose, (3.9) 

'^r^\t=tc,r=rc 
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= ^ (1 + Zc)sm0, (3.10) 

where 6 describes the direction of the null geodesic in the cluster rest frame (see Fig. |2|). 
Note that Eq. ( p.^ ) is a condition to fix the degree of freedom for the affine transformation. 
It is easy to see that Eq. (3.4) is satisfied by these initial conditions. The redshift at each 



point on the null geodesic is defined by 

l + z = -t'. (3.11) 




At the intersection between the past light-cone of the cluster of galaxies located at 
z = Zc and the decoupling hypersurface, we stop integrating Eqs. (|3.lD-(^^) and label the 
redshift at this intersection as zd(^;-2c)- Here note that zd(0;Zc) = z^. Then, the CMB 
temperature Tci{6;Zc) observed at the cluster of galaxies located at z = Zc is obtained by 

Tciie; Zc) = T-^t-^^D (rc (9; Zc)) , (3.12) 
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where r = rj){9;Zc) is the radial coordinate of the intersection between the past hght-cone 
of this cluster and the decoupling hypersurface. If we know the temperature 7b (r) on 
the decoupling hypersurface, we find Tci{9;zc) by the above equation. Conversely, if we 
know the information about Tci{9; Zc) through the kSZ effect, we find 7b(r) for the domain 

rB{TT;Zc) <r< rD(0;Zc). 

3.1 Case of homogeneous decoupling hypersurface 



As mentioned in the second to the last paragraph in § |2.2| , we do not have enough informa- 
tion to determine the physical quantities on the decoupling hypersurface, 7b (r), r]i){r) and 
(XD{f) and hence the embedding of the decoupling hypersurface t = t-D{r). To determine 
them, we need to impose two additional assumptions on these quantities. Here, we assume 
Td^v) and r]D{r) are constant, i.e., T^^r) = and r]-£,{r) = r]^. Then ai){r) also becomes 
a constant determined by Eq. ( p. 25 ). Consequently, from Eq. ( p.24| ), we can determine the 



decoupling hypersurface t = t^lr). 

Since an LTB universe model is neither homogeneous nor isotropic for an observer 
at an off-center cluster of galaxies, not only a dipole component but also higher multi- 
pole components in Tci{9; Zc) may exist. However, we find from our numerical results that 
the anisotropy in Tci{9] Zc) is dominated by dipole, and thus we can fit Tci{9; Zc) with the 
Doppler shifted temperature distribution, 

Tdipoic(e; Zc) := ^^^?^^r,(0; z,), (3.13) 

C — Vc{Zc) COS 9 

where Vc{zc) is regarded as the effective drift velocity of the cluster at z = z^ relative to 
the CMB rest frame, and is determined by using the least-square method (see Fig. ^). 

The result is shown in Fig. ^ compared with the observed values for nine clusters 
reported in |34, 5|, 36 1. The observational data are listed in Table H. Obviously, the 



Table 1: Observed drift velocities for nine clusters. 



i{l - 9): name 


redshift(2;i) 


drift velocity [km/s](wi) 


c7+[km/s] 


0-. [km/s] 


1 


A1689|5] 


0.18 


+170 


+815 


-630 


2 


A2163|J 


.41 




0.20 


+490 


-hl370 


-880 


3 


A226l|5t] 


0.22 


-1575 


-M500 


-975 


4 


A2396|5|] 


0.23 


-hl900 


-F6225 


-2650 


5 


A1835|5|] 


0.25 


-175 


-M675 


-1275 


6: Zw 3146 


35 


1 


0.29 


-400 


+3700 


-1925 


7: RX J1347-1145||] 


0.45 


+U20 


-M170 


-1270 


8: CI 0016 + 16|5| 


0.55 


-4100 


+2650 


-1625 


9: MS 0451 


3[ 


>1 


0.55 


+490 


+1370 


-880 



effective drift velocity in our model is far larger than the observed values. This fact means 
that our LTB universe model whose Tj){r), r]D{r) and aD{r) are constant cannot explain 
the observational data. 
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cos(e) 

Figure 3: The temperature anisotropy seen at a cluster of galaxies, Tc\{0; Zc), and the dipole fit to 
it, 7dipoie(^; Zc), for several values of Zc- 



Here we have assumed that the ratios of two of any physical quantities (the energy 
densities of baryon, dark matter, photons, temperature of photon) are constant on the 
decoupling hypersurface. Since, by virtue of the assumption ( p.lOD , our background LTB 
universe model is very close to a homogeneous Einstein-de Sitter universe at the decoupling 
epoch, these assumptions imply that inhomogeneities are small and adiabatic at that time. 
Therefore the above result strongly indicates that anti-Copernican LTB cosmological mod- 
els are virtually excluded by the observational data within the context of the conventional 
adiabatic perturbation scenario for the structure formation. 



In |]17|, homogeneity at the last scattering surface was also assumed and their 



results are qualitatively consistent with the result obtained in this section. 
3.2 Case of inhomogeneous decoupling hypersurface 

We now consider the possibility of an inhomogeneous LSS and look for an r-dependent 
decoupling hypersurface t = t]:){r) and the physical quantities T^^r), r]Y){r) and a-£){r) 
so that the observational results of the kSZ effect can be explained by our LTB universe 
model. 



First, for simplicity, we assume r/D(^) = i]* (see [25] about inhomogeneity of r/o). Then, 
to fix the decoupling epoch, we consider a comoving fiducial cluster of galaxies at redshift 
Zc = Zf. We assume that the anisotropy of the CMB temperature observed in the rest 
frame of the fiducial cluster is purely dipolar, i.e., 

Tci{9;zi) = M O;zi) = ^^^(l + Zf)To, (3.14) 

c — V{ cos c — V{ cos 8 
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Figure 4: Effective drift velocity at each cluster redshift for a homogeneous temperature distribu- 
tion on the decoupling hypersurface. 



where Uf is the effective drift velocity of the fiducial cluster of galaxies relative to the CMB 
rest frame. Then, the decoupling hypersurface, t = t-£,{r) is obtained from the above CMB 
temperature Tc\{0\Zi) in the following manner: We numerically integrate null geodesic 
equations (|3.1| )- (|3.3D pastward from this fiducial cluster of galaxies to every direction. At 
each redshift z on the null geodesic to the direction specified by 0, we define T and a by 

r= i±^rei(0;zf), (3.15) 
1 + 2f 

Po \T J 

where t = tg{z; 6) and r = rg{z; 6) represent the trajectory of this null geodesic (see Fig. |5|). 



During the numerical integration, we check at each redshift z whether Eq. ( 2.25 ) is 



satisfied by T-q = T , a-Q = a and t^d = If Eq. ( |2.25 ) is satisfied, we stop integrating 



Eqs. ( ^.lD -( |3^) . At this moment, we have z = z-£,{9;z{) and rg{z;9) = rD{6;z{). This de- 
termines T and a as functions of r: T = TD(rD(6'; Z{)) and a = 00(^0(6'; Z{)), respectively. 
The schematic figure for the trajectory r = rj:,{6; Z{) is shown in Fig. ^. We find from this 
figure that 

TBiTT] Zf) < r^ie; Zf) < rD(0; Zf) = n^{z,) for < ^ < vr. (3.17) 

The above inequality means that the functional forms of 7b (r) and aY){r) are determined 
only for the range rD(vr;Zf) < r < nc(z*). However, this information is enough for the 
calculation of the temperature anisotropies from other clusters at redshift Zc < Z{, because 
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the past light-cone emanated from this cluster intersects with the decoupling hypersurface 
inside the region rD(vr; Zf) < r < ric{z^) as shown in Fig. ^. 



LSS for Fiducial Cluster 



LSS for a Cluster 



LSS for the observer 




Figure 5: Schematic figure of LSSs on the decoupling hypersurface for the fiducial cluster, a cluster 
at a redshift Zc < Z[ and the observer eX z = (us). 



Applying the same procedure as in the previous section, we obtain the CMB temper- 
ature Tci{6; Zc) observed at a cluster of galaxies located at z = Zc < Z{. Results are shown 
in Fig. ^, where we set Zf = 0.6. The effective velocities fc(^c) for various values of V{ are 
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Figure 6: The effective velocities as functions of the cluster redshift z^ are shown for various values 
of Vf. The fiducial cluster is at Zf = 0.6. 
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shown as functions of the cluster redshift z^. The effective velocity of a cluster is evaluated 
from the temperature anisotropy seen from the cluster using the least-square method. It 
is clear that the result has much better consistency with the observational data than the 



model with an adiabatic situation in S3.1 



In order to compare our result with that obtained in [17|, we have calculated the same 
log- likelihood as theirs, 

-21n£ = ^i^i^#l±^. (3.18) 

As shown in the = 1600km/s, even though our model has an extremely large 

void, it turns out to be consistent with the observational data in the l-a level. 

However, it is too early to conclude from the above result that the LTB-type anti- 
Copernican cosmology is consistent with observation. There is a very important point 
to be reminded. In Fig. ^ we depict the temperature and density distributions on the 
decoupling hypersurface, T-£,{r) and pT){r) := p{tTi{r),r). We can see from Fig. 8 that there 
is a factor of two difference in the matter density poif) on the decoupling hypersurface 
over the range of r of our interest, while the temperature T£){r) is nearly uniform. Thus the 
uniform temperature hypersurfaces are far different from uniform density hyper surf aces, 
implying the existence of an isocurvature perturbation with the amplitude of order unity. 
In particular, on a p =const. hypersurface, the distribution of the radiation is highly 
inhomogeneous. Therefore it would be wrong to assume the evolution of the universe to 
be the same as in a homogeneous and isotropic universe in the radiation dominated era. 

The above result leads to a significant difficulty in predicting almost all observables 
related to the structure formation. To begin with, it implies that an LTB model in which 
the photon energy density is neglected ceases to be a good approximation to the universe at 
the decoupling epoch. Therefore we first have to obtain spherically symmetric background 
universe models with radiation and dust, in which the radial profiles of radiation and 
dust are diff'erent. Then we need to perform cosmological perturbation analyses of these 
models in order to make theoretical predictions. This is an issue well beyond the scope 
of this paper. We just mention that there have been some advancements in this direction 
recently [|^, |2|]. 



4. Summary and Discussion 

In this paper, we have discussed the kSZ effect in an LTB universe model whose distance- 
redshift relation observed at the center agrees with the concordance ACDM model at red- 
shift z <2. Since the LTB universe model is spherically symmetric, it is natural to assume 
that physical quantities at the decoupling epoch are functions of the radial coordinate. We 
have assumed the decoupling between baryons and photon occurs instantaneously, or in 
other words, it occurs on a single spacelike hypersurface. To determine this decoupling 
hypersurface we have adopted Gamow's criterion, H = T, for simplicity, where H is the 
Hubble parameter and F is the rate of collisions between photons and electrons. How- 
ever, this criterion is not sufficient to determine the decoupling hypersurface. Instead, we 
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Figure 7: Log-likclihoods for (crpv,t;sys) = (400,0), (400,750), (1600,0) and (1600,750). 



have found that we need two more conditions for the physical quantities (the energy den- 
sities of dark matter, baryons and photons, the temperature of photon) on the decouphng 
hypersurface. 

In § |3.1| , we have assumed the adiabatic perturbation scenario, namely we have fixed the 
ratios between two of any physical quantities (e.g., baryon-to-photon ratio) to be constant 
on the decoupling hypersurface. Then, it turned out that the CMB dipole anisotropy 
in each cluster rest frame is much larger than the observational data deduced from the 
kSZ effect. Therefore we conclude that our LTB universe model in the context of the 
conventional adiabatic perturbation scenario is ruled out. 
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Figure 8: r-dependence of Tj^ir) and PBir)- 



Then in §3.2, we have introduced radial non-adiabatic (isocurvature) inhomogeneities 
in the non-relativistic matter on the decouphng hypersurface. We assumed that the baryon- 
to-photon ratio is constant on the decoupling hypersurface. Then, we have shown that 
appropriate spatial distributions of isocurvature inhomogeneities can improve the consis- 
tency with observational results. It should be noted that we have introduced only a radial 
inhomogeneity. Hence it does not contribute directly to the CMB temperature anisotropy 
observed at the center. Actually, we may have a much larger value of the likelihood than 
the result in |17|. Thus, at present, the LTB universe model has not yet been ruled out 
completely by the observation of the kSZ effect. However, we must be cautioned by the fact 
that the necessary magnitude of isocurvature inhomogeneities is so large that our simple 
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LTB universe model that neglects the photon energy density ceases to be a good approx- 
imation to the universe at the decoupling. Therefore our result should be taken not as 
evidence but as an indication that the presence of large isocurvature perturbations may 
save the LTB universe model. 

Finally, we mention that the possibility to consider the off-center observer in LTB 
universe models. In [11|, Alnes and Amarzguioui claimed that the observer has to be 
located within a radius of 15Mpc from the center for the induced dipole to be less than that 
observed by the COBE satellite. Kodama et al. also reported similar constraints on several 
LTB models However, as shown in §^]^, it is possible to make an off-center observer 
in an LTB universe see a vanishing dipole anisotropy if we turn on radial isocurvature 
inhomogeneities at the decoupling epoch. It may be also interesting to investigate the 
relation of the kSZ effect in LTB universe models to the large scale coherent bulk flow 
reported by Kashlinsky et al. |38, 39, 4C]. 
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